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An analytical expression of the strain distribution due to lattice mismatch is obtained in an 
infinite isotropic elastic medium (a matrix) with a three-dimensional polyhedron-shaped inclusion 
(a quantum dot). The expression was obtained utilizing the analogy between electrostatic and elastic 
theory problems. The main idea lies in similarity of behavior of point charge electric field and the 
strain field induced by point inclusion in the matrix. This opens a way to simplify the structure of 
the expression for the strain tensor. In the solution, the strain distribution consists of contributions 
related to faces and edges of the inclusion. A contribution of each face is proportional to the solid 
angle at which the face is seen from the point where the strain is calculated. A contribution of an 
edge is proportional to the electrostatic potential which would be induced by this edge if it is charged 
with a constant linear charge density. The solution is valid for the case of inclusion having the same 
elastic constants as the matrix. Our method can be applied also to the case of semi-infinite matrix 
with a free surface. Three particular cases of the general solution are considered — for inclusions of 
pyramidal, truncated pyramidal, and "hut-cluster" shape. In these cases considerable simplification 
was achieved in comparison with previously published solutions. A generalization of the obtained 
solution to the case of anisotropic media is discussed. 

PACS numbers: 68.65.Hb, 46.25.-y 



I. INTRODUCTION 

Self-assembled quantum dots are three-dimensional in- 
clusions of one material in another one (a matrix). Usu- 
ally there is a lattice mismatch between materials of 
an inclusion and a matrix. The lattice mismatch gives 
rise to a built-in inhomogeneous elastic strain which in 
turn produces significant changes in the electronic band 
structure.— 1^ Therefore knowledge of the strain distri- 
bution is of crucial importance for electronic structure 
studies. Nearly all papers concerning electronic struc- 
ture calculations of quantum dots start with evaluation 
of elastic strain. Especially important is the strain dis- 
tribution for type-II quantum dots where the confining 
potential for one type of carriers is mainly due to the 
strain inhomogeneity."^ 

There are a lot of theoretical works on the strain 
distribution in quantum dot structures (for a review, 
see Refs. |^i5|). In addition to numerical calculations 
(using finite difference ^lii^ finite elementj^iiS valence 
force field^'^ d^d^d'^ and molecular dynamics^'* methods), 
some analytical techniques have been proposed. Most of 
them are based on the usage of Green's functions, ei- 
ther in the real spac o^^d^d"^ or in the reciprocal space. 
Some authors break the inclusion into infinitely small 
"bricks"— or into infinitely thin cuboidsi^ and then ap- 
ply the superposition principle. For ellipsoidal inclusions, 
Eshelby's approach^" has proved to be effective. Also 
a number of results obtained in thermoelasticity theory 
may be applied to lattice-mismatched heterostructures, 
as pointed out in Ref. [2l|. 

Different methods have their own merits and restric- 
tions. In our opinion, an ideal solution of the elastic in- 



clusion problem has to be analytical, to be expressed in 
terms of elementary functions and written in closed form, 
to be applicable to a broad range of inclusion shapes, 
and to take into account elastic anisotropy and atomistic 
corrections. Analytical closed-form solutions have been 
found for few cases of inclusion shapes: an ellipsoid,— 
a cuboid,^- a pyramidfi5.'22 g^^^j variety of quantum- 
wire-like structures 1^ Nozaki and Taya^'' have presented 
a general solution for an arbitrary polyhedron, but it is 
extremely complicated. All these solutions imply elas- 
tic isotropy and (except the case of ellipsoidal inclusion) 
equal elastic constants of the two media. 

The aim of our paper is to develop a novel approach to 
constructing the solutions for the general case of a poly- 
hedral inclusion, and to propose a new insight into the 
structure of a solution. We stress that the solution should 
have a clear physical or geometrical meaning. Without 
having a clear structure of a solution, it is hardly possible 
to develop its generalization to anisotropic media and/or 
to inclusions with clastic constants different from ones of 
the matrix. 

This paper considers the following problem. There is 
an infinite elastically isotropic medium (a matrix) with a 
finite polyhedron-shaped inclusion. The crystal lattice of 
the inclusion matches the lattice of the matrix without 
any defects. Elastic moduli of the inclusion are assumed 
to be equal to ones of the matrix, but the matrix and the 
inclusion have different lattice constants. This produces 
an elastic strain in both the inclusion and the matrix, and 
the task is to determine the strain tensor as a function of 
coordinates, eap{v). We neglect atomistic and nonlincar- 
ity effects, assuming that the lattice mismatch is small, 
and lattice constants are small in comparison with the 



2 



inclusion size. 

It is important to note that the strain distribution pro- 
duced by an inclusion in a semi-infinite matrix may eas- 
ily be calculated, provided that the corresponding strain 
field in an infinite matrix is known. We will discuss it 
in Section WUl 

The rest of the paper is organized in the following way. 
In Section [TTl a new approach to evaluation the strain 
distribution, based on an analogy between electrostatic 
and elastic problems, is described. The solution for an 
arbitrary polyhedron-shaped inclusion in an infinite ma- 
trix is presented and discussed in Section IIIII Then, in 
Section IIVI this solution is applied to pyramidal, trun- 
cated pyramidal, and "hut-cluster" inclusions. Section FVl 
shows the possibility of generalization of our method to 
anisotropic media. Section IVII contains the summary of 
the paper. The Appendix is devoted to evaluation of solid 
angles that is important for calculation of the strain. 



II. ELECTROSTATIC ANALOGY 
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FIG. 1: A sketch of modification of the "charge distribution" 
by taking derivatives of the potential: (a) transformation of 
a volume charge into a surface charge, (b) transformation of 
a surface charge into a linear charge and a dipole layer. 



The starting point of our investigation is a well-known 
analogy between the elastic inclusion problem and the 
electrostatic problem (Poisson equation).-- Namely, the 
displacement vector u(r) induced by the inclusion is pro- 
portional to the electric field F(r) that would appear if 
the inclusion were uniformly charged: 



u(r) 



£o(l- 



47r(l - V 



■F(r) 



£o(l- 



47r(l - v) 



r — r' 



'|3 



dr', (1) 



where Eq is the lattice mismatch (eq — (oinciusion ~ 
amatrix)/amatrix, a being the lattice constant), v is the 
Poisson ratio, V denotes volume of the inclusion. For 
simplicity, in our auxiliary electrostatic problem we take 
the charge density and the dielectric constant equal to 
unity. Zero displacements correspond to positions of 
atoms exactly in sites of the ideal lattice of the matrix. 
Strain tensor is defined as 



dxR dxa 



£0(50/3 x(r), (2) 



where Xa is a-th component of the position vector r. Sap 
is the Kroneker delta, and x(r) is equal to 1 inside the 
inclusion and to outside it. 

Introducing an electrostatic potential 



dr' 



we can express the field F(r) as 

F(r) = -V^(r). 
Combination of Eqs. ([T]), ([2]) and (|4|) produces 



-A 



dxadxp 



£0 5al3X{^), 



(3) 



(4) 



(5) 



where A = eo(l -I- i>)/Att{1 - v). 

Now our aim is to evaluate second derivatives of the 
potential </5(r). For this purpose, we introduce three ad- 
ditional functions: 

1) Ti{Y) — an electrostatic potential of the uniformly 
charged (with unit surface density) z-th face of the inclu- 
sion surface; 

2) $fc(r) — an electrostatic potential of the uniformly 
charged (with unit linear density) fc-th edge of the inclu- 
sion surface; 

3) Q.i (r) — an electrostatic potential of the dipole layer 
uniformly spread over the i-th face with surface density 
of dipole moment equal to n* — the outward normal to 
the face. 

The potential $(r) of an uniformly charged edge is 
expressed as an integral / dl / \y — y'\, where dl is a linear 
element of the edge, and r' is a position vector of this 
linear element. Evaluation of this integral gives: 



$(r) = log 



ri+r2 



ri + r2 



(6) 



where ri and r2 are distances from the point r to ends of 
the edge, and L is the edge length. 

The potential r^(r) of a fiat uniform dipole layer is 
known2& to be equal to the solid angle at which this layer 
is seen from the point r, taken with positive sign if the 
positively charged side of the layer is seen from r, and 
with negative sign otherwise. Thus below we will refer to 
the quantity rii(r) as to a solid angle subtended by the 
i-th face from the point r. 

In order to find second derivatives d'^tp(Y) j dxadxp, we 
note that the first derivative can be expressed as a sum 
over faces of the inclusion surface: 



9r 



(7) 



(faces) 



(see Fig. la). Indeed, this derivative can be rewritten as 



dr 



d_ 
dr 



X(r') 



r - r' 



dr' 



9x(r') dr' 



d{r') |r-r'| 



(8) 



The derivative 9x(r')/9r' plays the role of a "charge den- 
sity" in Eq. ([S]). It vanishes everywhere except the sur- 
face of the inclusion. Near the i-th face of the surface, 
x(r) is equal to — 0(n*(r — r^)), where is any point 
of this face, 9 is the Heaviside function; consequently 
dx{^)ldr = — n*^(n*(r — r^)), that corresponds to the 
"surface charge density" — n* at the i-th face. Therefore, 
a contribution of the i-th face to d(p{r)/dr is equal to 
— n*^i, according to Eq. ([7]). 

The next step is finding of derivatives of !Fi{r). Let, 
for simplicity, the i-th face lie in the plane XY, and its 
outward normal vector be directed along the axis Z. 
To find derivatives dTi/dx and dTijdy, one can follow 
the same line of argumentation as at deriving Eq. ([7]). 
The result is: 



dx 



dy 



(9) 



where k runs over edges surrounding the i-th face; b*^ is 
a unit vector which is parallel to the i-th face, directed 
out of this face, and perpendicular to the fc-th edge (see 
Fig. lb). The last derivative, d!Fi/dz, transforms an uni- 
formly charhed z-th face into a dipole layer with surface 
density of dipole momentum equal to — n*, as shown in 
Fig. lb. Consequently, 



m 

dz 



(10) 



Eqs. ^ and pH)) can be written together in a vector 
form, which is independent on orientation of a face with 
respect to co-ordinate axes: 



dr 



-n'a(r)-Eb'=$fe(r) 



(11) 



Using Eqs. and ([TT]), one can express second deriva- 
tives of the potential (p{r) via solid angles fli{r) and po- 
tentials of charged edges $fc(r): 



aV(r) 

dxa dx 



E 

i 

(faces' 

E 

k 

(edges 



n>*5 n,{r) 

{nfbf+n^^bf)^u{v), (12) 



where for each edge k there are four unit vectors n 



fei 



,fe2 



related to the two faces intersecting at 
this edge. With given normals n'^^ and n'^^, the vectors 
b*^^ and b*"'^ can be found in the following way: 



where l'^ is a unit vector directed along edge k. From two 
possible directions of 1*^ , one should choose the one going 
clockwise with respect to face kl and, correspondingly, 
counter-clockwise with respect to fc2, when seeing from 
the outside of the inclusion. 



III. GENERAL SOLUTION AND ITS 
PROPERTIES 

A. The general solution 

In the previous Section, it was shown that the strain 
tensor eaj3{r) can be expressed (via Eq. (O) in terms of 
second derivatives of some auxiliary "electrostatic poten- 
tial" 93(r). In turn, these second derivatives break down 
into contributions of all faces and edges of the inclusion 
surface (Eq. (fT^ ). Combining equations (O and (|12p . 
we obtain the following expression for the strain tensor: 



e„/3(r) = - A E ">^)f^*(i-) 



- A E 7a/3*fe(r)-£o'5a,3X(r), (14) 

k 

(edges) 

where i runs over faces of the inclusion surface, and k 
runs over its edges. 

In Eq. p4)) . ili(r) is a solid angle subtended by the j-th 
face from the point r (positive if the outer side of the face 
is seen from the point r, and negative otherwise); <I'fc(r) 
is the electrostatic potential of an uniformly charged fc-th 
edge (with unit linear charge density) at the point r; x(r) 
is equal to 1 inside the inclusion and to outside it; eo is 
the relative lattice mismatch between the inclusion and 
the matrix; the constant A is equal to eQ{\-\-v)/An{l — v), 
where v is the Poisson ratio; 6ap is the Kroneker delta; n* 
is a normal unit vector to the i-th face, directed outside 
the inclusion; and a constant tensor 7^^ is equal to 



7a/3 



fcl>fel 



fe2.fe2 
T-a 0/3 ■ 



(15) 



In Eq. (Uni), n*^! and n*^^ are normal unit vectors (directed 
outside the inclusion) to the two faces which intersect at 
the fc-th edge; b*^^ is a unit vector perpendicular to the 
k-th edge and to n*^^, and directed out of the fcl-th face; 
analogously, b*^^ is a unit vector perpendicular to the fc- 
th edge and to n*^^, directed out of the fc2-th face (see 

Eq. inD). 

The tensor 7^^ is symmetrical, and it can be written 
in an equivalent form 



7a/3 



= K4-i?^i3^) sin0. 



where 
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FIG. 2: Cuboidal inclusion. Numbers in circles refer to faces, 
italic numbers refer to edges. 



and 9 is the internal dihedral angle between faces fcl and 
k2. 

There is a simple expression ^ for potentials <&fc(r) 
contributing into Eq. (|14p . Some closed- form expressions 
for solid angles fli{r) are presented in the Appendix. 

Equation (|14p is the main result of the present paper. 
It gives a closed-form analytical expression for strain dis- 
tribution in and around a polyhedral inclusion buried 
into infinite isotropic elastic medium. 

With a known strain tensor, one can easily obtain the 
stress tensor a^is via Hooke's law)^ 



E 



1 



£Q/3(r) 



(16) 



where E is the Young modulus. 

On the basis of Eq. p^ . the program code is written 
that can easily calculate the strain distribution produced 
by a lattice-mismatched inclusion in an infinite, isotropic 
matrix. Inclusion shape can be an arbitrary polyhedron. 
This program, named "easystrain" , is freely available at 
pittp : / / easystrain . narod . ru 



B. Cuboidal inclusion 



C. Hydrostatic strain 

Now we consider some simple properties of the solution 
([T^ . These properties can be regarded as tests of validity 
of the solution. 

First, let us calculate the hydrostatic component of 
strain (that is, the trace £aa(r) of the strain tensor). 
Taking into account that n^n^ — (n*)^ = 1, = 
j^fei|-,fei _^ j^fc2|-,fc2 _ g^j^^ _ 3^ readily get from 



Eq. Ill] 



«(r) 



3£ox(i-) 



The sum of solid angles, ^^i(r), vanishes for any point 
r outside the inclusion. Indeed, all faces can be di- 
vided into two groups with regard to the point r: 1) the 
ones whose outer sides are seen from the point r, 2) the 
ones whose inner sides are seen from r. Net solid an- 
gles subtended by the two groups are the same, but they 
contribute to the sum with opposite signs and 

therefore cancel each other. 

If the point r is inside the inclusion, all the faces belong 
to the second group and the solid angle subtended by 
them together are the full solid angle, 4tt. So, J^i ^ j(r) = 
—An. Combining both cases (r outside and inside the 
inclusion) we get 



^a(r) = -4^x(r), 



(17) 



and consequently 

£o„(r) = (47rA - 3eo)x(r) = -2eo^-^x{r)- (18) 

1 — 1^ 

We have come to the well-known result that the hydro- 
static strain is zero outside the inclusion and is constant 
inside iti^i 



If the inclusion has the form of cuboid with faces per- 
pendicular to the direction of the axes x, y and z (Fig. [2]), 
then separate components of Eq. are simplified to 

e^^(r) = -A(0i(r) + f^2(r))-£ox(r), 
e,y{r) = -A($i(r)-$2(r) + $3(r)-$4(r)), 

and all the other components of strain tensor have a sim- 
ilar form. So, the diagonal components Exx-, Syy and Ezz 
depend only on solid angles fii, whereas off-diagonal com- 
ponents Exy, Exz and Eyz depend only on edge contribu- 
tions $fc. This is in agreement with results of Downes, 
Faux and O'Reilly siS These authors pointed out that, in 
the case of cuboidal inclusion, solid angles subtended by 
faces contribute into the stress tensor (and hence into the 
strain tensor too) . Our paper generalizes this observation 
to the case of any polyhedral inclusion. 



D. Strain discontinuities at faces 

Then, it is easy to examine the behavior of the strain 
at the inclusion surface, starting from Eq. When 
the point r, moving from outside toward inside, crosses 
a face of the inclusion, the strain changes stepwise. The 
discontinuity of the strain, Ae^^ = e"|"^° — e°'^*''"^°, can 
be written as 

AEai3 = -A^n>^Ar2i - A^7^^A$fe 

i k 

There Afi^ and A$fe denote discontinuities of fi^ and '^k- 
In fact, edge contributions $fc have no discontinuities at 
the face, and only one of fi^ has a discontinuity — namely, 
the Q,i related to the face under consideration. For this 
face, Afii = — 47r, because the solid angle reaches 
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FIG. 3: Quantum wire inclusion. Point A is the point where 
the strain is to be obtained. 6i and axe the angle and 
distance contributing to Eq. (|21|l . 



27r at the face, and the value Vli changes the sign from 
positive to negative. So, 

with a unit vector n perpendicular to the face. It is 
convenient to consider the strain tensor with respect to 
the coordinate axes rj, ( connected to the face: axes f 
and 77 parallel to the face, and the axis C perpendicular 
to it. Therefore 71^ = n,, = 0, rt^; = 1, and Asap takes 
the following form: 

Ae^^ = Ae,„ = -eq, Ae^c = ^ZT^' ^^^^ 

off-diagonal components Ae^^, Ae^^^, Ae^ii^ arc zero. Ac- 
cording to Hooke's law (fTB|) . discontinuities of the stress 
tensor, Actq/j, are 

Aa^5 = Aa,^ = -j^, ^^CC = (20) 

and again off-diagonal components are zero. 

Equations p9|l and ([20|) are in accordance to the 
boundary conditions at the interface: Ae^^ = Ae^^ = 
— £0, Ae^,j — (continuity of displacement field u(r)), 
Aa^Q ~ Aa^j^ ~ Aa^^ — (balance of elastic forces at 
the interface). 

E. Quantum-wire inclusion 

Next, we consider the strain distribution in a quantum- 
wire-like inclusion and its surrounding (Fig. [3]) . Such an 
inclusion is a prism, the top and bottom of which go to 
infinity. For simplicity, let all side faces and edges be 
parallel to the axis z. So the strain is independent on z. 



To obtain the strain distribution, one may start from 
Eq. (|14p for a prism of finite height, and then go to the 
limit of infinitely large vertical dimension. In this limit, 
contributions of base and bottom faces, as well as of edges 
adjoining to these faces, vanish. For each side face, the 
solid angle reduces to a doubled plane angle 9i (Fig.[3|) 
subtended by the cross-section of this face by a plane 
parallel to axes x, y: fli = 2di. Edge contributions $fc 
reduce to simple logarithmic expressions: 

<I>i: = — 2 logdi + const, 

where dk is the distance to the fc-th edge (Fig. [3]). The 
constants in these expressions are infinitely large, but 
they cancel each other being substituted into Eq. (fT4|) . 

As a result, we come to the following expression for the 
strain in a quantum- wire-like inclusion: 

Safsir) = - 2A^7i>^6ii 

i 

+ 2A^7^^1ogdfc -£o<5a/3X- (21) 

k 

There the indices i and k run over all side faces and edges, 
correspondingly; 9i is the plane angle subtended by the 
cross section of i-th face by the plane passing through the 
point r parallel to the axes x and y (positive if the outer 
side of the face is seen from the point r, and negative 
otherwise); dk is the distance from the point r to k-th 
edge. All the rest notations are the same as in Eq. (|14p . 
As the XZ-, yz- and zz-components of the tensors n^an^p 
and 7^^ are zero, the corresponding components of strain 
tensor are independent on 9i and d^: 

Equation (j2ip is an equivalent, but more simple and 
compact, form of the solution obtained by Faux, Downes 
and O'Reilly,^ 

F. Semi-infinite matrix 

Finally we discuss the strain distribution in a semi- 
infinite matrix. Davies^S, proposed a method of reducing 
the elastic inclusion problem in a semi-infinite matrix 
to the corresponding problem in an infinite matrix. For 
convenience, we reproduce there the results of Davies's 
workjS£ 

Let an inclusion be buried in a semi-infinite matrix that 
fill a half-space z > Zg, or z < z^, with a free surface in 
the plane z — Zg. Isotropic linear elasticity is assumed, 
and elastic moduli of the matrix and the inclusion are 
the same. To calculate the strain distribution £0/3 (r) in 
this system, one can previously find an analogous strain 
distribution Sap (r) in a system consisting of the same in- 
clusion in an infinite matrix. It can be found by Eq. (I14p . 
for example. Then, components of eQ;/3(r) are expressed 
via components of £a^(r), ea/3(r2) and their derivatives 
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dsai3/dz{r2), where the point r2 is a "mirror image" of 
the point r with respect to the surface: 



Y2 



{x,y,2zs - z), 



Syy{V 

£yz{r 



£xx{r) + (3 - 4iy)e^.^{r2) + 2{z - z,) 



Br 

^^xx 

dz 



(r2), 



de 



Syyiv) + (3 - ^v)eyy{V2) + 2{z - Z s) ^ [V 2) , 

rip 

e^r) - (1 - 4i.)e,,{r2) + 2{z - ^s)^(r2), 

de 

s,y{r) + (3 - 4i/)e,y(r2) + 2(z - z,)-^{r2), 
£xz{r) - exz{r2) - 2(z - z^)— ^(r2), 



Syz{r) - £yz{r2) - 2(z - Zg) 



dz 
dey 
dz 



(r2). 



If the inclusion is a polyhedron, Eq. (|14p provides an 
analytical expression for the strain tensor e^^. Therefore 
its derivative dsap /dz can be evaluated analytically as a 
combination of derivatives of solid angles 17^ and values 
It is important to note that derivatives of fli can be 
expressed via derivatives of 



dn^jr) 

dXa 



dxp 



(22) 



where summation is over all the edges adjoining to the 
i-th face; and b*^* is the one of unit vectors b*^^, b'"'^, 
which is perpendicular to n'. Eq. may be useful 
since analytical expressions for solid angles are rather 
complicated in comparison with the expression ([6]) for 
values 



IV. APPLICATION TO PYRAMIDAL AND 
HUT-CLUSTER INCLUSIONS 

Among all polyhedrons, the three ones appear most 
often as geometrical models of quantum dots. These are 
square-based pyramid (Fig. 3^) , truncated square-based 
pyramid (Fig. dJ^), and so-called "hut-cluster" (Fig. O;). 
In this Section, we apply the general expression ^1)1 to 
the specific cases of pyramidal and hut-cluster inclusions. 
The case of truncated pyramid does not demand a spe- 
cial consideration, because it is easy to obtain solution for 
truncated pyramid, provided that the solution for pyra- 
mid has yet been obtained (see below). 



A. Pyramid 

With the numbering scheme of Fig. 0^, we get the fol- 
lowing expressions for tensors n^n^ and 7^^: 

n>o = {0,0, 1,0,0,0}, 
nafip = {s^ 0, 0, -sc, 0} , 




FIG. 4: Inclusions of most common shapes: (a) a pyramid 
with square base; (b) a truncated pyramid with square base; 
(c) a "hut-cluster". Numbers in circles refer to faces (the 
base has the number 0), italic numbers refer to edges, and 
bold numbers — to vertices. 
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{0, s", c2, 0, 0, -sc] , 
{s^, 0, c^, 0, sc, 0}, 
{0, s\ c\ 0, 0, sc}, 
s X {c, 0, — c, 0, — s, 0} , 
s X {0, c, — c, 0, 0, — s} , 
s X {c, 0, — c, 0, s, 0} , 
s X {0, c, — c, 0, 0, s} , 



'1- 


hc2 


'l+c2 


'1- 


h c2 



{— c^, — c^, 2c^, 1, sc, sc} , 
{-c^, -c^, 2c^, -1, -sc, sc] 
{— c^, — c^, 2c^, 1, — sc, — sc} 

-c^ 2c^ 



sc, 

1, sc, — sc} , 



where s = sini?, c = cosi?, and i? is a dihedral angle 
between the pyramid base and any of its side face. The 
tensor components are listed in braces in the following 
order: xx, yy, zz, xy, xz, yz. 

It is worth to note the following property of the set of 
tensors 7^0: 



k 



0, 



where is the length of the fc-th edge. This property 
comes from a requirement that all terms proportional to 
r^^ in Eq. \1A\ must cancel each other at r — s- 00. It may 
serve as a useful test of correctness of the results. 
These values of nj^n^ and 7^^, together with Eq. 
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give the expression for the strain distribution in a pyra- midal inclusion and its surrounding: 



J 



^yy 

£xy 
Syz 



s^Ai^i + n^) - scA($i + $3) + 7f=f *5-8 - Sox, 

^5-8 " £oX, 



s2 

-s^A{n2 + Qi) ~ SCA($2 + $4) 
-AOo - c^Arii_4 + scA$i_4 - 

^^^^ .(*5-*6 + *7-$8), 

sc' 



sc^A 



A(l]i - ria) + .s2A($i - $3) - ^'^^ 
2 



^$5-8 -£oX, 



^^($5-$6-*7 + $8), 



SCA(1]2 - f^4) + s"A($2 - $4) - 7=f (*5 + $6 - *7 - 3'8)- 

I 



(23) 



There we use a shorthand notation: rii_4 — fli + 0.2 + 
r^s + U4, and so on. Note that, using Eq. (fT7)l. one can 
simphfy the expression for Szz to the foUowing one: 



Ezz = -s^Afio + scA$i_4 



2se^A 



5-8 



eo(c^i±^-l)x. 



(24) 



Eq. IjASp in the Appendix provides analytical expressions 
for solid angles r2o---f^4 in the pyramid. 



B. Truncated pyramid 



To get the solution for the truncated pyramid, ^q™"'^'', 
the easiest way is to start from the solution for a 

pyramid, e^J^^^ and apply the superposition principle. 
The full pyramid, OABCD, consists of a truncated 
one, ABCDAiBiCiDi, and a smaller one, OAiBiCiDi 
(Fig. |3|d). According to the superposition principle. 



-(pyr) I 

'-al3 



£ 



(trunc) 

Q/3 



(r) 



_{small) 

'a/3 



(r), 



where e^J'"'' ^i*™"'^^ ^^'^ ^i™^"^ reier to figures OABCD, 
ABCDAiBiCiDi and OAiBiCiDi, correspondingly. 
Then, it is well known that, in the framework of the con- 
tinual elasticity theory, similar inclusions produce similar 
strain fields. As pyramids OABCD and OAiBiCiDi are 
similar, there is a relation between e^J'^ ^^'^ ^ap^^^^^- 



.(pyr) , 

'a/3 ' 



_{small) 

-a/3 



(ro + Ar), 



where tq is a position vector of the apex O, A is a trunca- 
tion parameter (a ratio of sizes of the two pyramids, see 



Fig. Sh). So e^j'fl""'^'' can be expressed in terms of e)^ 



(pyr) . 



e 



(trunc) 
a/3 ' 



r) 



-(pyr; 

-a/3 



(r) 



.(pyr) 

-a/3 



A 



+ ro). (25) 



This solution was compared numerically with the so- 
lution published in Ref. |22| . We reproduce all strain pro- 
files presented in that paper, except the component e^z 



in Fig. 11 of Ref. |22, where the absolute value coincides 
with our results, but the sign was opposite. We believe 
that the sign of Sxz in Ref. [2^ is erroneous, because it 
leads to an incorrect behavior of the strain at r — s- cx). 
Indeed, the multipole expansion, being applied to Eq. (fT|). 
gives for large r 



-a/3 



(r + i-c) = 



47r(l — v) 



+ 0(r-5). 



where is a volume of the inclusion, and Tc is a position 
vector of its center of mass. For a pyramid, Xc = Uc = 0- 
This expression shows that, at fixed positive x and y, 
£xz must be negative when z +00 and positive when 
z ^ — 00. This predicted behavior of Sxz disagrees with 
Fig. 11 of Ref. 122, but agrees with our calculations. 



C. Hut-cluster 

Finally, we consider the hut-cluster. The hut-cluster 
is a figure that consists of the base (a parallelogram) 
and four side faces. Slope angles of all the side faces 
are the same. Therefore orientations of faces and edges 
of the hut-cluster are the same as of pyramid, except 
the top (9th) edge. So, the solution for the hut-cluster 
is very similar to the one for the pyramid. The only 
difference is the addition of the contribution of 9th edge. 
Of course, the values of solid angles Qo---^i and of edge 
contributions $i...<l>8 in the hut-cluster are not the same 
as in the pyramid. Analytical expressions for solid angles 
in the hut-cluster are given by Eq. (|A6p in the Appendix. 

To get the solution for the hut-cluster from Eq. ([23]) , it 
is sufficient to add the term 2scA$g to Syy, and to add the 
term — 2scA$g to Ezz- This demonstrates the flexibility 
of the general solution (fTi)) . This is a property that is 
not inherent in previous particular solutionsJ^iiSii^S 

An analytical solution for a hut-cluster was first ob- 
tained by GlasAS as a special case of a more general an- 
swer for a truncated pyramid with rectangular bottom 
and top faces. Our method provides a much more simple 
solution. 
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D. Strain profiles along axes of symmetry 

If the point r lies at the four-fold axis of symmetry 
of the pyramid, expressions are simplified greatly, 
because all values r2i(r)...ri4(r) are the same, values 
$i(r)...$4(r) are the same, and $5(r)...$8(r) are also 
the same. Moreover, it is sufficient to evaluate only the 
2;2:-component of the strain, because non-diagonal com- 
ponents Exy, Sxz, Syz are zero, and other diagonal compo- 
nents £xx and Eyy can be expressed via Szz using Eq. (fT8)l : 



-yy 



2 i^aa 



1 - 2i/ 



I 



1 



V 2 



(26) 



It is convenient to use Eq. ([2l| for evaluation the com- 
ponent Ezz- Let us put the origin of coordinate system 
to the center of the pyramid base. So, coordinates of 
a point at the axis of symmetry are (0,0, z). Let h be 
a pyramid height, a be a length of the pyramid base, 
I = 1/0^/2 -I- be a length of a side edge, and r be a 
distance between the point (0, 0, z) and any vertex of the 
pyramid base: 



According to Eq. (jA5p . the solid angle Hq is equal to 



Then, 



fig = —4 arctan 
$1 = $2 = $3 = ^'4 = log 



4zr 



2r + a 



2r - a 

\z-h\+r + l 



\z ~ h\ + r - I' 
Substituting all these quantities into Eq. ((24)) . we get 



8ahA f 2h 2r + a 

e^^(0,0,z) = 2 I ^^2 — arctan —+ log 

a-^ + 4/i^ V a Azr 2r — a 



a \z — h\ -'r r -\- 1 
1 °^\z~h\+r-l 



Ax- 



(27) 



Here we expressed quantities s and c in terms of a, h, 
and I; A = eo(l -I- i^)/47r(l - i/); x = 1 if 2; S (0; /i) and 
X — Q otherwise; the constant A is the coefficient at x in 
Eq. jH: 



a^(l+i^) _ 1 
(a^+4/i^)(l-i^) ^ 



A^eo{c''j^-l) = eo 

For the truncated pyramid, 
Eq. ^ give 



Eq. together with 



ezz{0,0,z 
+ log 



8ahA 
{2r+a){2f~a) 



2h 

— arctan 

a Azr 



2h 

— arctan — rr 
a Azr 



a Xr + r+{l-X)l 
1 °^Af + r-(l-A)Z 



(2r-a)(2f+a) C"'' Xf + r-il-\)l J ^^^^ 
In Eq. dig), X = 1 if ^ e (0; (1 - A)fe), an d x = oth- 
erwise; z = {z — h) / \ -\- h; f = 1/0^72"+^ is a distance 
between the point (0,0, z) and any vertex of the base; 
I = y/a'^/2 + h?. 

In a similar manner, one can get the strain profile along 
the axis of symmetry of the hut-cluster. As there is only 
two-fold axis in the hut-cluster, the components £xx and 
eyy are no longer the same. Therefore we cannot use 



Eq. (|26p to extract £xx and Syy from Ezz- Instead, we 
should find Szz and £xx independently, and then extract 
Eyy by means of Eq. : 



1 - 2v 
^-2eo- X 



1 



We chose the center of the cluster base as an origin of 
the coordinate system. Let h he a cluster height, a and 
b be the smaller and the bigger edge lengths of the base, 
correspondingly. Then, values of Ezz and Exx at the axis 
of symmetry are 



„ , 8ahA f2h ab 

,{0,0, z)^ —arctan- 

a'^+Ah'^ \ a Azri 



1 {2r,+a){2r, + b) r, + r, + l _l 2r,Hh~a) \ . 

2 ^(2ri-a)(2ri-&) I ^ r^ + r^-l 2 ^ 2r^-{b-a) } ^ ^ ' 



ExxiO,0,z) 



8ahA 
i^+Ah^ 



2h a^z — abh 1 2ri+a 

— arctan [^(^_„)^4^(2_;,)]^^^(^2^.4^2)^j. ~2 2^^ 



— log 

21 r^+n-l 



r 



-eoX, (30) 



where ri (z) = 
the point (0, 

r,{z) ^ y/{b^a)yA+{z 



y/a^/A + W/A + z^ is a distance from 
z) to the first vertex of the hut-cluster; 

ft.) 2 is a distance from the 



point (0, 0, z) to the fifth vertex; / — y/a^j2 + h^ is a 
length of each side edge; x = 1 if G (0; ft) and x = 



otherwise. The last log term in Eq. is a contribu- 
tion of the 9th edge, and the arctangent term in Eq. ((30|) 
comes from the first and the third faces. All the rest 
terms are similar to that of Eq. (|27p . 

As an illustration, in Fig. [5] we plotted profiles of the 
strain component, Ezz, calculated by Equations 
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FIG. 5; The strain component, e^^, plotted along the z axis 
for a pyramid (Eq. (|27|l . solid line), a truncated pyramid 
(Eq. ([28j, dash line), and a hut-cluster (Eq. lf29|l . dash-dot 
line) . 

for some particular cases of pyramidal, truncated pyra- 
midal, and hut-cluster inclusions. Parameters of the 
structures are chosen to be the same as ones in Ref. [l^: 
eo = 0.067; v = 0.24; a = 155 A; /i = 55 A for the 
pyramid and the hut-cluster; h = 110 A and A — 0.5 for 
the truncated pyramid; b = 2a for the hut-cluster. Pre- 
sented curves for the pyramid and the truncated pyra- 
mid are identical to ones of Ref. [13 (curves D and B in 
Fig. 5, correspondingly), that confirms the correctness of 
our formulas. 

A detailed discussion of these profiles is beyond the 
scope of the present paper. We only note that Szz di- 
verges logarithmically at 2 = 55 A in the pyramid and 
the hut-cluster. This divergence is a common feature of 
a strain distribution in a vicinity of a vertex or an edge 
of any polyhedral inclusion. 

V. ELASTIC ANISOTROPY 

The above consideration was based on an assumption 
of elastically isotropic inclusion and matrix. For applica- 
tions to semiconductor heterostructures, this assumption 
may be a source of considerable error. For example, the 
Young modulus of silicon in (111) direction is 1.44 times 
greater than that in (100) direction. Therefore, taking 
the elastic anisotropy into account is an actual problem. 

In this Section, we argue that our method can be ex- 
panded to anisotropic media. We start from expression of 
strain tensor via Green's tensor Gaf3 by Faux and Pear- 
son3^: 

Eapir) = £0 / Ga/3(r ^ r')dr' , 
Jv 

where V is the inclusion volume, and Eq is lattice mis- 
match. These authors found a series expansion for 
Green's tensor, assuming cubic anisotropy: 

Ga/5 = GS + AG«-f A2gS + ..., (31) 



where expansion coefficient A = (Gn — G12 — 
2C44) / {C12 + 2G44) is a measure of anisotropy (A ss — i 
for typical semiconductors), Gn, G12 and G44 are elastic 
moduli. Each term of this expansion can be presented 
as a combination of partial derivatives of expressions like 
1/r, x^/r, etc. For example, isotropic term is 

^(0), s ^ _£o 3Gi2 + 2G44 1 
"1^^^' 47r G12 + 2G44 dxc^dxpr' 

xx-component of first-order correction, gS, is a linear 
combination of the following terms: 

I X x^ qa p 

dx"^ r' dx^ r' dx'^ r ' dx^dy"^ r ' dx'^dz^ r 

As a result, strain tensor £Q^(r) expresses as a combina- 
tion of derivatives 

where a, 6, c, d,e, f ~ 0, 1, 2, ... with a constraint a + b + 
c = d + e + f + 2. Each term in the expansion ((3T|) is 
a sum of a finite number of derivatives (j32p taken with 
proper constant coefficients. 

According to our method, the integrals in Eq. ([32]) can 
be regarded as electrostatic potentials induced by a non- 
uniformly charged inclusion. Taking the derivatives in 
Eq. ([32|) , one proceeds from "volume charge" to "surface 
dipoles" on faces of the inclusion surface, and to "lin- 
ear charges" and "multipoles" on its edges. Thus, our 
method allows to split the strain tensor into contribu- 
tions of faces and edges of inclusion surface (assuming 
that inclusion shape is a polyhedron): 

e„Mr)=^A«(r)+^i?S(r), (33) 

i k 

where indices i and k run over all faces and edges, corre- 
spondingly. 

. . (k) 

Explicit formulas for contributions A^j^ and B^^ are 
beyond the scope of the present paper and are the subject 
of a separate publication. There we only note that each 
face contribution A^*^(r) is proportional to a sohd angle 
fli(r) with a coefficient depending on the orientation of 
this face and on elastic constants. Edge contributions 
b'^P (r) can be expressed in a closed form for each term 
of series expansion (PT|) . 

VI. CONCLUSIONS 

In summary, we propose a new, more simple and flexi- 
ble expression for strain field in and around an inclusion 
buried in an infinite or semi-infinite isotropic medium. 
This expression was also implemented as a computer 
programi^ We show that the strain field can be pre- 
sented as a sum of contributions of the faces and edges. 
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FIG. 6; Solid angles: (a) expressed by Eq. (|A2|l : (b) expressed 
by Eq. (|X4|). 



This is the main point of our method; it gives a possi- 
bility to construct expressions for strain distribution in 
inclusions of complicated shapes. The general solution is 
applied to important particular cases of pyramidal and 
hut-cluster inclusions. Our solution for the pyramid re- 
produces previous solutions, but in a simpler and intu- 
itively understandable form. We believe that it paves the 
way for further simplifications and generalizations of the 
solution, for example, to the case of anisotropic elasticity. 
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APPENDIX A: 



EVALUATION OF SOLID 
ANGLES 



Here we present some explicit formulas expressing solid 
angles as functions of coordinates. 

A solid angle r2(r), that a surface S subtends at a point 
r, may be defined as an integral over the surface: 



(Al) 



where dS is a surface element, n is a unit vector directed 
normally to this surface element, and r' is a position 
vector of the surface element. 



This integral is easily evaluated if the surface S' is a 
rectangle. For simplicity, let this rectangle lie in the plane 
z = 0, and its edges be oriented along the axes x and 
y (Fig. E^). Let xi and X2 be ^-coordinates of edges 
directed along the axis y {xi < X2)', yi and y2 be y- 
coordinates of the rest two edges of the rectangle (yi < 
2/2). Then the integral (jAip is expressed as follows: 



Vt'""'^{r;xi,X2,yi,y2) ^ dx' 



dy' z 



{{x-x'y^ + {y-y'y + z^f/^' 



xi VI 



or 



n''"'\r;xi,x2,yi,y2) 
arctan (^-^^)iy-y^) 

z rii 

arctan (---^Ky-v^ + 

Z 1-21 



_ arctan ii^ii^ia^ 

z ri2 
ix-x2){y-y2) 

2 r22 



arctan ■ 



(A2) 

Here rii...r22 are distances from the point r to the corners 
of the rectangle: 



'•2l=V {x~X2Y + {y-y\Y+z'^ , r22 = ^ (x-X2Y + {y-y2Y+z^. 

It is assumed that values of arctangents fall into the range 
(-§>+§)• 

To find a solid angle subtended by a triangle (Fig.lHb), 
one can use the relation by Oosterom and Strackeei^Si^ 



triangle 



tan- 



(a,b,c) 



[a X b] c 



abc+ (ab)c+ (ac)6+ (bc)a 



(A3) 



Vectors a, b, c join the point O, at which the solid 
angle is subtended, to the vertices A, B, C of the triangle. 
To satisfy the condition that the solid angle is positive 
if it is looked at from outside, and is negative otherwise, 
one should choose the proper order of following of the 
vertices A, B, C. Namely, the closed contour ABCA 
should follow in a clockwise direction, seeing from outside 
the inclusion (as shown by dashed arrows in Fig.[6lD). So 
the triple scalar product [a x b] c is positive (negative) if 
the outer (inner) side of the triangular face is seen from 
the point O. 

Care must be taken while resolving Eq. (jA3p with re- 
spect to fl. For simplicity, we will refer to the right side of 
Eq. (jASp as to A. The sign of fl is the same as the sign of 
the product [ax b] c, but may differ from the sign of A. So 
we cannot "naively" resolve Eq. (jASp as = 2 arctan A. 
Instead, we should write down = 2 (arctan A mod tt) if 
[a X b] c > 0, and = 2(arctan A mod vr) — 2tt otherwise. 
Joining together both cases, we obtain 



^triangle (j^^^^j.) 

= 2 ( arctan 



[a X b] c 



abc+ (ab)c-j- {a.cjb+ (bc)a 
- 27r(l-6i([ax b] c)) , 



mod TT 



(A4) 



where d{x) is the Heaviside function (1 for positive x, 
for negative x). Note that the easiest way to implement 
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Eq. (|A4p in a computer program is to use C math library 
function atan2: — 2* atan2(P, Q), where P and Q are 
numerator and denominator of the right part of Eq. (|A3p . 

With Eq. (|A4p . one can calculate a solid angle sub- 
tended by an arbitrary polygon, breaking this polygon 
down into triangles. 

Now let us apply the expressions (|A2|) and (jA4p to 
faces of the pyramid and the hut-cluster considered in 
Section HVl First, we choose a reference frame with the 
origin at the center of the pyramid base, the axes x and 
y along edges of the base, and the axis z directed to- 
ward the apex of the pyramid. So, position vectors of 
the vertices (see Fig. |4^) are 

ri - (a/2, a/2,0), = (-a/2, a/2, 0), 

ra = (-a/2, -a/2, 0), = (a/2, -a/2, 0), = (0, 0, h), 

where a is a base edge length, h is a height of the pyra- 
mid. Note that a dihedral angle — arctan(2/i/a). Solid 
angles contributing into Eq. are 







(r;- 


a a 
2 ' 2 ' 


a a\ 
2' 2'' 






^triang 


^'=(ri 


- r, r4 


- r,r5 


-r), 


f^2(r) 


^triang 


^'=(r2 


- r, ri 


- r,r5 


-r), 




_ ^triang 




— I", I"2 


- r,r5 






_ ^triang 


'^(r4 


- r,r3 


- r,r5 


-r). 



The minus sign at 17''°'^* in Eq. (|A5|) reflects the fact that 
the pyramid base is directed downwards. 

For the hut-cluster, position vectors of vertices are 

n = (6/2, a/2,0), r2 = (-6/2, a/2,0). 



ra = (-6/2, -a/2,0), = (6/2, -a/2, 0), 
1-5 = {{b - a)/2, 0, h), re = (-(6 - a)/2, 0, h). 

Here a and 6 are the smaller and the bigger edge 
lengths of the base. (Again, the dihedral angle is ■& — 
arctan(2ft,/a).) We introduce also two points Oi and O2 
where side edges cross (see Fig. |4}d): 

roi = (0, -(6 - a)/2, hb/a), ro, = (0, (6 - a)/2, hb/a). 

The solid angle ^2 of the trapezoidal 2nd face is the 
difference of two solid angles subtended by triangles Oil2 
and Oi56. The angle ^4 is evaluated in the same way. As 
a result, solid angles subtended by faces of the hut-cluster 
are 



1^0 (r) 


= f7''°"*(r: ^ 


b a 
2 ' 2 


2 




J^i(r) 


= f^trianglc^j.^ _ 


r,r4 - 


r,r5 - 


-r), 


f^2(r) 


^ ^triangle ^j.^ _ 


r,ri - 


i",roi 


-r) 






r,r5 - 


r,roi 


-r), (A6) 




= ^^'"^"^'^(ra - 


r,r2 - 


r, re - 








r,r3 - 




-r) 




- ^2*"^"si°(r5 - 


i",r6 - 


r,ro2 


-r). 



These expressions can be substituted into a modified ver- 
sion of Eq. , as described in Section HVl 
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